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The mass of the W boson is one of the least precisely measured parameters of the electroweak interaction. Confronted 
with other measurements of standard model parameters it can test the internal consistency of the Standard Model and 
can constrain the possible mass of the standard model Higgs boson. The CDF collaboration has published the current 
QQ , single most precise measurement of the W boson mass using 200 pb of CDF Run II pp data, and an improved 

C " 3 ' measurement with 2.4 fb~^ of CDF Run II data is underway. Both measurements are described in detail in these 

' proceedings. 

(N : 

00 1. INTRODUCTION 

m . 

__, The standard model of particle physics (SM) has been very successful in describing electroweak interactions. The 
^ mass of the W boson is one of the parameters of the electroweak theory, which at tree level is fully determined by the 
^ relations to the other electroweak parameters. The loop corrections to the W boson mass calculation are dominated 
by the loop containing the b and t quarks and could receive a contribution from a loop with the yet-unconfirmed Higgs 
boson, while further corrections can come from loops with possible supersymmetric particles. Precision measurements 
of the W boson mass allow us to constrain the size of the possible loop corrections, giving us valuable information 
T-H on the range of possible Higgs masses. The proceedings present one of the most precise single measurements of the 
^ W boson mass. 

\o . 

P 2. THE W BOSON MASS MEASUREMENT AT CDF 

a^ ■ 
o : 

OO i The Tevatron experiments CDF and D0 have long been involved in the measurement of the W boson mass, first 

^P. I pushing the precision below 100 MeV/c^, and later published the combined average of the Run I W boson mass with 

' uncertainty of 59 MeV/c^ [1]. In 2001 the upgraded CDF and D0 detectors began taking data and the size of the 
Run I dataset has long been surpassed. The CDF detector uses its tracking and calorimetry capabilities to measure 

^ , the W boson mass both inW ev and W fiiy decays. 



2.1 . Measurement Strategy 

The W bosons in pp collisions are mainly produced in quark anti-quark annihilation. In the environment of hadronic 
machines, only W ev and W ^ ^iv decays can be reconstructed with the purity needed for the W boson mass 
measurement. The neutrino is not detected in the CDF detector and its momentum cannot be fully kinematically 
constrained, since the momentum of the incoming partons producing the W boson is not known, and the momentum 
along the beam-line cannot be accurately measured. The momentum conservation of the measured quantities can 
instead be applied in the transverse plane, since the partons have negligible transverse momenta. The transverse 
momentum of the neutrino (p^) is thus inferred from IfT^ the missing transverse momentum required to achieve 
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Figure 1: (a) The fractional momentum correction Ap/p as a function of the mean inverse momentum of the muons from 
J/i/j T(15') — > s-nd ^ ~^ decays. The full line is the best linear fit to the points, and the dotted lines 

mark the total (statistical and systematic) uncertainty of the fit. (b) The E/pc distribution used for the calibration of the 
energy measurement. 



the momentum balance in the transverse plane. The so-called transverse mass is constructed out of the kinematic 
quantities in the transverse plane: 

niT = ^J 2plj. • • [1 - cos{(f>i - 0^)] . (1) 

The dependence of rriT on mw is determined using Monte Carlo (MC) templates, which are created for different 
input TOvF ■ The template best corresponding to the data distributions is extracted by a binned maximum- likelihood 
fit. The W boson mass is also obtained from the templates constructed for lepton transverse momentum plj^ and the 
transverse momentum of the neutrino, I^t- The final mw value is a weighted average of the three determinations. 
While the three variables are highly correlated, they have different sensitivities to boosts and systematic effects, 
hence their combination improves the mw determination. 

For the desired accuracy of the mw measurement the tot, Pti templates need to be generated using a 

far greater number of events than would be feasible with full detector simulation. A fast parameterized detector 
simulation is developed for this purpose and is described in detail in Ref. [2] . The fast event generation and detector 
simulation enables reliable estimates of systematic uncertainties since we are able to produce several simulation 
samples where the relevant parameters are varied within their uncertainty. 

The main quantities in the event measured by the detector are the momentum of the charged lepton and the 
hadronic recoil, the total transverse momentum of the remaining particles, recoiling (in the transverse plane) from 
the W boson. The measurement of the recoil is necessary to reconstruct the missing transverse momentum I^t- 
Recoiling particles come from gluon radiation of the incoming quarks, and the underlying event processes. The fast 
event generation and detector simulation have to be calibrated to predict these kinematic variables to a part in 10^. 



2.2. Momentum and Energy Measurement Calibrations 

The muon momentum is determined using the Central Outer Tracker (COT), a gas- filled cylindrical drift chamber 
placed in a f .4 T magnetic field. The internal alignment of the COT is performed using high-pT cosmic rays. The 
fast simulation of the momentum measurement takes into account the energy loss and scattering of particles in the 
material, COT hit resolutions and the track-fitting algorithm. The absolute momentum measurement scale and hit 
resolutions as well as the number of simulated radiation lengths of traversed material are tuned on data. 

The energy of the electron is obtained from the measurement in the electromagnetic calorimeter (EM), while the 
angular information of the electron track is obtained with the COT. The simulation of the energy measurement 




Figure 2: (a) The comparison of the simulated magnitude of the recoil with data, for W ^ ev candidates, (b) The comparison 
of the simulated distribution of the recoil component in the direction of the charged lepton with data, for W —> fiv candidates. 



includes the nonlinearity of the calorimeter response, the energy dependence of the energy measurement and the 
shower leakage out of the EM calorimeter. 

The calibration of the measured momenta and energies is one of the most important aspects of the W boson mass 
measurement, since scaling the measured momenta directly shifts the fitted kinematic distributions. The absolute 
momentum scale is set by requiring that the measured values of the precisely known masses of the Z, J/ip and T(IS') 
resonances in the decays to a pair of muons correspond to their PDG values [3]. Figure 1(a) shows that the fractional 
momentum scale corrections obtained from different invariant mass fits are consistent with each other. Combining 
the calibration using all three particle masses enables momentum determination with an accuracy of 0.021%. 

Taking advantage of the high J/^p — > fJ-'^fi~ and T(15') fJ-^fJ-^ statistics, the measured momentum can be cali- 
brated more accurately than the energy measured in the calorimeter. The momentum calibration is hence transferred 
to the energy measurement using the E/pc distribution (Figure 1(b)). We adjust the absolute energy scale on the 
peak of the E/pc distribution: the tuning is done in several transverse energy bins to account for non- linear calorime- 
ter response. The number of radiation lengths of material used in the energy loss simulation is tuned to match the 
radiative tail of the E/pc distribution. 

The calibration yields a, Z e+e~ mass measurement of Mz — 91190 ± 67stat MeV/c^, in very good agreement 
with the world average (91187.6 ±2.1 MeV/c^ [3]). The calorimeter calibration is further improved by combining the 
tunes of the E/pc and the Z — » e^e~ invariant mass distributions. Similar to the scale, the resolution of the energy 
measurement is tuned on the widths of the E/pc and the Z — > e+e^ mass distributions. 



2.3. Hadronic Recoil Calibration 

Simulation of the reconstructed transverse momentum of the neutrino, i.e. the missing transverse momentum I^t, 
requires an accurate simulation of the hadronic recoil, the total transverse momentum of all particles except that of 
the two prompt leptons from the W boson decay. It is measured from all the energy depositions in the calorimeter, 
except those associated with the primary charged lepton. These depositions are caused by particles from initial state 
QCD radiation, the hadronization of the pp remnants not partaking in the boson creation, multiple collisions from 
the same bunch crossing and from initial and final-state photons. The contribution of different components to the 
magnitude and spread of the recoil distributions is parameterized and the parameterization tuned on the Z £^£~ 
decays and minimum bias data (for the underlying event). The comparison of the data and the simulated magnitude 
of the recoil (Figure 2(a)) and the component of the recoil in the lepton direction (Figure 2(b)) shows good 
agreement. The mw measurement is specially sensitive to U|| as it directly affects the calculation of mx- 



2.4. Generation of W boson production and decay 



The dedicated parameterized Monte Carlo simulation has to model both initial state gluon radiation, which is 
responsible for a non-zero W boson pt, and the QED radiation of the final state lepton(s). The W boson pt 
spectrum is modelled with RESBOS [4] and its parameterization of the non-perturbative low pT region is tuned on 
the dilepton px distribution in the Z i^fr decays. The QED corrections of onc-photon emission arc simulated 
using WGRAD [5]. The photon energies are increased by 10% (with an absolute uncertainty of 5%) to account for 
additional energy loss due to two-photon radiation which is not included in the WGRAD calculation. 

Systematic uncertainties due to the limited knowledge of the non-perturbative QCD parameters determining the 
W Pt spectrum and the simulation of the radiation of more than one photon were estimated (see Table II). 

Parton distribution fimctions (PDFs) arc used to calculate the kinematics of the W boson decay products. The 
CTEQ6M [6] PDFs and their error sets are used for the generation of templates and the PDF uncertainty estima- 
tion (Table II). 

2.5. Bacltgrounds 

Due to imperfect detector measurement and coverage, several background processes can mimic a W ^ £1/ decay, 
thus passing the W event selection. We estimate the amounts and shapes of the dominant background contributions 
and add them to the simulation templates when performing the mw fits. 

The main backgrounds in the W ^ ev channel are Z e+e" and W ^ tv events, and events where a hadronic jet 
is mis-reconstructed as a lepton ('multi-jet background'). We model Z e+e^ and W ^ tv events using Pythia [7] 
and the full GEANT simulation of the CDF detector. The multi-jet background normalizations are found by fitting 
the low distribution, where this background dominates. 

The main backgrounds in the W ^ /j.v channel are Z /^"""M" and W ^ tu events, cosmic rays, decays of low 
momentum pions or kaons into muons in the COT. and hadronic jets, penetrating into the muon chambers or decaying 
to muons. The Z ^ fi'^IJ'~ and W ^ tv events are modelled using Pythia and the full detector simulation. The 
contribution of low momentum meson decays is found by fitting the high tail of the track fit distribution where 
this background is large. 

3. W BOSON IVIASS FITS AND THE RESULT 

Finally, the W boson mass is determined separately for the electron and muon channels, in the fits to the three 
kinematic distributions (tot, Pt -^t)- The results of the fits are summarized in Table I. The two transverse 
mass fits are shown in Figure 3. The breakdown of the systematic uncertainties for the example of the tot fit is 
given in Table II. 

The six W boson mass fits are combined, taking into account the correlations, to give mw = 80413 ± 34stat ± 

34syst MeV/c^. With a total uncertainty of 48 MeV/c^, this measurement is the most precise single measurement of 
mw to date. This result increases the world average W boson mass to mw = 80398 ± 25 MeV/c^ [8], reducing its 
total uncertainty by 15%. 

A global electroweak fit with the new mw average together with the current best measurement of the top quark 
mass, mt = 172.6 ± 1.4 GeV/c^ [9], constrains the possible Higgs boson mass to mn — 86^27 GeV/c^ [8]. The 
corresponding 95% confidence level upper limit on the Higgs boson mass is 160 GcV/c^ [8]. Figure 4 shows the 
comparison of the mw and mt measurements with their constraint on the Higgs boson mass. The mw and mt units 
of the plot refiect how the uncertainty on both masses propagates to the uncertainty on the predicted mn- To achieve 
the same level of constraint from W boson and top-quark mass measurements requires their uncertainties to satisfy 
Amw ~ 0.006 ■ Am(. It can be seen from Figure 4 that currently the uncertainty of the mw measurement is the 
limiting factor, thus motivating further work. 




The total uncertainty of the pubhshed CDF measurement is smaUer than it was projected for 200 pb~^ from the 
Run I CDF and D0 measurements, mainly due to the inclusion of the T(IS') and J/ip resonances in the momentum 
calibration and the usage of the E/pc distribution for the energy measurement calibration. Further improvements 
in the detector model and the production and decay model (for example QED radiative corrections) are likely to 
reduce other systematic uncertainties as well. The CDF collaboration therefore predicts that the measurement of 
the W boson mass with a precision better than that of the current world average (25 MeV/c^) is possible with the 
CDF data already in hand. 



4. FIRST LOOK AT THE 2.4 fb-^ CDF RUN-II DATA 

The CDF collaboration has now collected well over ten times more data than were used in the published analysis [2] . 
Most of the systematic uncertainties in the published measurement (Table II) are limited by the statistics of the 
calibration samples. The sensitivity of a precision measurement to much higher instantaneous luminosities and the 
fact that the data from a much longer period of data-taking have to be calibrated, prevent a straight-forward inclusion 
of more data. 

The CDF collaboration has selected a data sample corresponding to approximately 2.4 fb~^ that is most suitable 



Distribution 


mw (MeV/c^) xVdof 


mriW eu) 
P^t{W eu) 
p'i{W eu) 


80493 ± 48stat ± 39,yst 86/48 
80451 ± 58stat ± 45syst 63/62 
80473 ± 57stat ± 54syst 63/62 




80349 ± 54,tat ± 27syst 59/48 
80321 ± 66stat ± 40syst 72/62 
80396 ± 66stat ± 46syst 44/62 



Table I: Fit results with total systematic and statistical uncertainties from all three fit distributions used to extract mw- 



for the improved mw measurement and have aheady performed several cahbrations with the data (due to different 
triggering requirements the W fiu and Z M"*"/^" decay samples correspond to approximately 2.3 fb"-'^, as stated 
on the plots). The following preliminary figures show that the increased luminosity and range of the data samples 
have neither significantly deteriorated the quality of the collected data nor our ability to measure the quantities of 
interest. 

4.1. Instantaneous Luminosity 

Figure 5(a) shows the comparison of the instantaneous luminosity distributions for the published analysis and the 
preliminary 2.4 fb~^ data. Higher average instantaneous luminosities and a much wider distribution of luminosity 
values affect the modelling of multiple pp interactions and backgrounds. The recoil distributions are the ones most 
severely altered and the recoil simulation now includes the instantaneous luminosity in its parameterization. The 
dependence of the hadronic recoil on the instantaneous luminosity can be seen in the distributions of 'SEt, the total 
recoil energy deposit in the calorimeter. Figure 5(b) shows the comparison of the data and the simulated distribution 
of T,Et for events with instantaneous luminosity smaller and greater than 70 x 10^° s~^ cm~^ (this value roughly 
corresponds to the average instantaneous luminosity of the whole sample) for Z — > e+e~ events. The agreement shows 
that the parameterizations are able to capture the dependence of the total recoil energy deposition on luminosity. 

4.2. Statistical scaling of uncertainties 

The degradation of the detector performance and the changes introduced to the detector over time are systematic 
effects that can spoil the scaling of measurement uncertainties with statistics. Special care has been taken with the 
data calibrations to keep the quality of the new data comparable to the one in the published analysis. 
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39 


27 


26 


Total Uncertainty 


62 


60 


26 



Table II: Systematic and total mw uncertainties for the transverse mass fits inW ev and W ij,v channels. The correlated 
part of each systematic uncertainty between the two channels is given in the last column. All values in MeV/c^. 
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Figure 4: The constraint on the Higgs boson mass from direct niw and mt measurements at 68% confidence level (CL). The 
band shows the SM relationship of the two masses as a function of itih which can be read from the abscissa. 
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Figure 5: (a) the comparison o the instantaneous luminosity distributions for the published analysis [2] and the preliminary 
2.4 fb~^ data, (b) The distribution of 'SEt for events at instantaneous luminosity greater and smaller than the average 
instantaneous luminosity of the preliminary 2.4 fb~^ dataset. The plots are for Z — » e"'"e~ events. 



Figures 6 and 7(a) show the invariant mass distributions used for the calibration of the momentum measurement. 
The fit for the momentum scale in J/ip — > IJ.~ decays is done in bins of average curvature, the example plot shown, 
Figure 6(a) is for the bin corresponding to the muons with the highest reconstructed momenta. Figures 7(b) and 8(b) 
show the kinematic variables used for the calibration of the energy measurement. The invariant mass of the candidate 
electrons from the Z e+e^ decays, using only the information from the COT, is shown in Figure 8(a). Due to 
the significant energy loss suffered by the electrons as they pass through the detector material this distribution has 
worse Z boson mass resolution than the invariant mass constructed from calorimeter information. Good agreement 
between data and simulation gives confidence in the simulation of the passage of electrons through the detector. 
Figure 9 shows the mx distributions for data and simulation, used in the myy fit- 

The statistical uncertainty of the relevant quantity is obtained from the fits. The statistical uncertainty of the 
mw fits of rriT and the mz fits of the invariant dilepton mass is obtained from a fit with the simulation blinded 



CDF II preliminary i.dr=2.3fb'^ CDF II preliminary Ldf=2.3fb"' 




(a) (b) 

Figure 6: The invariant dimuon mass for J/?/) — » fJ,'^ j-L~ (a) and T(15') f^^ decays (b). The plot for J/tp —> fJ,'^ is for 
muons with highest average momenta. The arrows indicate the limits of the fitting window. 
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Figure 7: The invariant dimuon mass for Z 
fitting window. 



fJ-^ H (a) and Z 



decays (b). The arrows indicate the limits of the 



by an unknown additive offset. In Table III the obtained uncertainties are compared to the expected statistical 
uncertainty obtained by scaling the published result [2] by the square root of the ratio of integrated luminosities. 
The fit uncertainties show that the resolutions of the distributions in the 2.4 fb~^ dataset are comparable to those 
of the data used in the published result. 



5. CONCLUSION 



We have published the most precise single measurement of the W boson mass [2]. The total uncertainty of the 
published CDF measurement is smaller than it was projected from previous measurements, and the improvements 
in the detector model and the production and decay model will further reduce the systematic uncertainties of future 
measurements. The CDF has set a goal of measuring the W boson mass with a precision better than that of the 
current world average (25 MeV/c^), using a dataset of approximately 2.4 fb^^. 

The preliminary plots using 2.4 fb^^ show that the expected statistical uncertainty of W mass fits scales with 
statistics, indicating that the potential degradation of the level of description due to the larger spread of instantaneous 
luminosities and time-dependent effects has not significantly affected the sensitivity of the fits. 
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Figure 8: (a) The invariant mass of the Z — > e^e electron pair, calculated from the momentum information measured in the 
COT. (b) The E/pc distribution of electrons from W ^ eu decays. The arrows indicate the limits of the fitting window. 
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Figure 9: The transverse mass distribution for W —> fiv (a) and W ^ ev decays (b). The arrows indicate the limits of the 
fitting window. 
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uncertainty 


published 


expected 


prelim, fit 


(all in MeV/c^) 


200 pb"^ 


2.3/2.4 fb"^ 


2.3/2.4 fb"^ 




43 


13 


12 


AmT'-iZ ^ e+e-) 


67 


20 


20 


Am|*"*-(Z^e+e-track) 


143 


42 


42 




54 


16 


16 


Aml\,^'-{W ev) 


48 


14 


15 




20 


6 


5 


Am^^'''''°*^* ''(mom. scale, J/V" — » /U"*"//"- one bin) 


20 


6 


6 


Am^*''*^'**''*'\mom. scale, T{1S) //"""/i") 


5 


1 


1 



Table III: Statistical scaling, all results in MoV/c^. The statistical uncertainty for the published result (200 pb~^), the expected 
statistical uncertainty obtained by scaling the published result by the square root of the ratio of integrated luminosities, and 
the statistical uncertainty obtained from the presented preliminary fit. 



